The  Determination  of  Wave  Group  Parameters  from  the  Resolution  of 
Incident  and  Reflected  Waves 


ABSTRACT 


Laboratory  tests  conducted  to  determine  the  effects  of  wave  groups  on  rubble 
mound  breakwater  stability  are  analyzed.  A  literature  review  outlining  the  controversy 
regarding  whether  or  not  wave  groups  effect  breakwater  stability  is  given.  Resolution  of 
this  controversy  is  proposed  using  two  independent  parameters  that  characterize  wave 
groups.  The  two  parameters  are  an  envelope  exceedance  coefficient  a  and  spectral 
shape  7.  The  incident  time  series  must  be  known  to  compute  a,  and  an  algorithm  is 
given  that  will  resolve  incident  and  reflected  wave  time  series  from  closely  spaced  wave 
gauges.  The  digital  to  analog  simulation  of  laboratory  waves  with  specified  wave  group 
characteristics  is  reviewed.  It  is  recommended  that  further  studies  be  conducted  to 
parameterize  the  envelope  exceedance  coefficient  a  with  a  non-biased  groupiness  factor. 
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INTRODUCTION 


The  design  of  rubble  mound  breakwaters  is  one  of  the  oldest  and  most  common 
applications  of  coastal  engineering  methodologies.  Design  of  rubble  mound  structures  is 
most  commonly  based  on  the  empirical  methods  presented  in  the  Shore  Protection 
Manual  (SPM  1984).  Due  to  some  failures  of  structures  designed  with  the  SPM 
methodology,  an  effort  to  improve  the  SPM  design  methodology  is  undertaken.  The 
design  methodology  presented  in  the  SPM  is  based  on  empirical  formulas  developed 
from  experiments  with  monochromatic  waves.  Most  of  the  more  recently  proposed 
design  criteria  are  based  on  laboratory  experiments  using  irregular  waves.  The  effects  of 
spectral  shape  and  wave  grouping  will  be  targeted.  As  a  proposed  improvement,  factors 
such  as  permeability,  wave  period,  and  storm  duration  may  have  secondary  effects  on 
rubble  mound  breakwater  stability. 

It  has  been  well-documented  that  waves  in  random  seas  tend  to  form  groups  not 
accounted  for  by  the  purely  Gaussian  random  wave  model.  Many  of  the  coastal 
engineering  problems  associated  with  wave  groups  have  been  documented  by  Medina 
and  Hudspeth  (1990).  In  the  last  two  decades,  controversy  has  arisen  as  to  whether 
spectra  composed  of  grouped  waves  cause  different  levels  of  damage  to  rubble  mound 
structures  than  spectra  of  ungrouped  waves.  Medina  et  al.  (1990)  have  found  that  most 
of  the  controversy  may  be  resolved  if  two  independent  spectral  related  parameters  are 
chosen  to  characterize  the  wave  groups  (Mase  and  Iwagaki  1986).  One  parameter 
characterizes  the  length  of  wave  groups  and  the  other  characterizes  the  magnitude  of  the 
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energy  flux  in  a  group  of  waves. 

Intermediate  scale  laboratory  studies  to  determine  if  wave  groups  affect  rubble 
mound  breakwater  stability  were  performed  in  the  2-dimensional  wave  channel  at  the  0. 
H.  Hindsdale  Wave  Research  Facility  (OHH-WRF)  at  Oregon  State  University.  In 
order  to  determine  the  effects  of  wave  groups  on  rubble  mound  structures,  sea  states 
with  identical  amplitude  spectra  but  different  wave  group  characteristics  were  simulated. 
The  resolution  of  incident  and  reflected  wave  time  series  from  the  random  wave 
simulation  is  required  in  order  to  evaluate  the  effects  of  wave  groups.  In  order  to 
quantify  the  two  wave  group  characteristics  proposed  by  Mase  and  Iwagaki  (1986),  the 
incident  time  series  must  be  known.  Medina  et  al.  (1994)  used  a  method  developed  by 
Kimura  (1985)  and  Fassardi  (1993)  that  extends  an  algorithm  proposed  by  Goda  and 
Suzuki  (1976)  to  resolve  the  incident  and  reflected  wave  time  series  from  closely  spaced 

wave  gauges. 
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CHAPTER  I:  Correlation  of  Wave  Groups  with  Armor  Damage 


Carstens,  et  al.  (1967),  first  observed  the  effect  of  wave  groups  on  breakwater 
damage  when  he  analyzed  two  time  series  with  significantly  different  groupiness,  but 
similar  spectra  and  wave  statistics.  The  time  series  with  the  more  grouped  waves,  taken 
from  a  field  site  in  the  Barents  Sea,  created  more  damage  than  the  time  series  generated 
from  the  theoretical  Neumann-spectrum.  Johnson,  et  al.  (1978),  as  a  result  of  their 
laboratory  studies  on  the  effects  of  wave  grouping  on  breakwater  stability,  concluded  that 
wave  grouping  was  an  essential  parameter  in  the  model  testing  of  rubble  mound 
breakwater  stability.  Wave  trains  w'ith  similar  spectra  and  wave  statistics,  but  different 
groupiness,  caused  significantly  different  levels  of  damage,  the  greater  damage  being 
caused  by  the  grouped  wave  trains.  Burcharth  (1979)  also  concluded  that  wave  grouping, 
or  the  succession  of  the  waves,  is  important  in  model  design. 

Since  the  energy  spectrum  does  not  contain  the  phases  of  the  individual  wave 
components  in  the  sea  surface  it  is  necessary  to  also  resolve  the  phases  of  the  individual 
wave  components  to  reproduce  an  accurate  succession  of  wave  heights  (Rye  1982).  In 
his  summary  of  breakwater  design,  Bruun  (1989)  recorded  that  the  resolution  of  the 
phase  spectrum  identified  apparent  order  where  none  is  assumed  in  the  Gaussian 
random  wave  model.  He  concluded  that  these  ordered  wave  groups  are  more  damaging 

than  other  waves. 

In  contrast,  van  der  Meer  (1988)  found  that  the  spectral  shape  and  grouping  of 
waves  had  little  or  no  influence  on  the  stability  of  breakwater  armor  layers  provided  that 


3 


the  average  spectral  period  was  used  to  simulate  the  sea  state  rather  than  the  peak 
period.  Hall  (1994),  in  a  study  on  bermed  breakwaters,  concluded  that  wave  groups  are 
only  important  during  the  stage  of  incipient  motion,  and  that  once  motion  is  initiated 
wave  groupiness  is  no  longer  important. 

Medina  et  al.  (1990)  found  that  much  of  the  controversy  could  be  resolved  if  the 
two  independent  parameters  of  wave  groups  proposed  by  Mase  and  Iwagaki  (1986)  were 
considered.  One  parameter  measures  the  run  length  or  the  number  of  high  waves  in  a 
time  series,  and  a  second  parameter  quantifies  the  magnitude  of  the  variation  of  wave 
energy.  In  their  laboratory  tests,  Medina  et  al.  (1990)  tested  four  different  wave  height 
time  series  proposed  by  Mase  and  Iwagaki  (1986)  shown  in  Figure  1.1,  while  keeping 
other  secondary  factors  that  affect  breakwater  stability  constant. 
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Figure  1. 1  -  Classification  of  wave  height  time  series  (Fassardi  1993) 
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Figure  1.1  shows  that  identical  spectra  may  have  different  magnitudes  of  wave  height 
variability,  which  may  cause  different  levels  of  damage.  Likewise,  different  spectra  may 
have  the  same  magnitude  of  wave  height  variability,  and  thus  cause  the  same  level  of 
damage.  Medina  et  al.  (1990)  simulated  truncated  Goda-JONSWAP  spectra  with 
different  run  lengths  and  different  magnitudes  of  variation  of  wave  energy  according  to 

Appendix  B. 

Medina  et  al.  (1990)  concluded  that  the  spectral  shape  and  run  length  alone  were 
not  sufficient  parameters  to  accurately  describe  wave  groupiness,  but  that  another 
parameter  that  measured  the  variation  of  wave  energy  more  completely  described  the 
wave  group.  Medina  et  al.  (1990)  correlated  an  envelope  exceedance  coefficient  or  that 
measures  the  variation  of  wave  energy  with  stability,  and  demonstrated  that  random 
wave  trains  with  the  same  spectral  shape  may  produce  different  levels  of  damage. 

Medina  et  al.  (1994)  found  no  significant  correlation  between  the  spectral  peakedness  7 
and  armor  damage,  but  concluded  that  the  envelope  exceedance  coefficient  a  was 
necessary'  to  evaluate  the  effects  of  random  waves  on  structures  in  both  physical  and 

prototypical  models. 
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CHAPTER  H:  Derivation  of  Envelope  Exceedance  Coefficient  or 


From  their  analysis  of  storm  waves,  Mase  and  Iwagaki  (1986)  recommended  using 
at  least  two  parameters  to  quantify  wave  groups,  one  parameter  that  measures  the 
magnitude  of  a  sequence  of  high  waves  in  a  time  series,  and  a  second  parameter  that 
measures  the  magnitude  of  the  variation  of  wave  energy.  Both  parameters  are  shown 
schematically  in  Figure  2.1,  where  A(x,t)  is  the  wave  envelope  function  from  the  Hilbert 
transform  defined  by  Medina  and  Hudspeth  (1990). 


Figure  2.1  -  Wave  group  parameters 
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The  peak  enhancement  factor  7  of  the  Goda-JONSWAP  spectrum  (Goda  1985) 
controls  the  sharpness  of  the  spectral  peak  and  was  selected  to  characterize  the 
magnitude  of  the  sequence  of  high  waves  in  the  time  series.  Rye  (1980)  related  the 
magnitude  of  the  sequence  of  high  waves,  run  length,  to  the  peak  enhancement  factor  7 
of  the  JONS  WAP  spectrum.  The  coefficient  a  was  selected  to  characterize  the 
magnitude  of  variation  of  wave  energy.  The  envelope  exceedance  coefficient  a  is  a 
dimensionless  measure  of  the  wave  energy  exceeding  that  associated  with  a  prescribed 
wave  height;  e.  g.,Hi/io. 

In  a  two  dimensional  wave  channel,  the  wave  height  function  H(x,t)  is  defined  by 

H(x,t)=  2A(x,t)  (2-^) 

If  At  =  a  constant  sampling  time  interval,  the  discrete  wave  height  function  at  a  fixed 
location  may  be  expressed  as  H(nAt).  Defining  a  normalized  measure  of  the  variation  of 
wave  height  as 

AH  =  1  (2.2) 

0  TT 


and  the  Heaviside  step  function  (Farlow  1982)  of  AH„  as 


U(AH„) 


0  AH^  <0 

1  ah”  >  0  ’ 

a 


(2.3) 


then  the  magnitude  of  the  variation  of  wave  energy  above  Hj/io  of  the  time  series  is  given 
as 
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(2.4) 


a'=  1  S  (AH/U(AH„) 
Nn=l 


where  N  is  the  number  of  discrete  data  points  in  the  time  series.  Normalizing  Eq.  (2.4) 
by  the  expected  value  of  a'  yields  the  envelope  exceedance  coefficient  a  defined  as 

(2.5) 

E[a'] 


where  E[q:']  is  defined  as 


E[a']  = 


H 

TT 

L  ^\no 


p(H)  dH 


(2.6) 


where  H  is  the  wave  height  function  at  a  fixed  location  H(t),  and  p(H)  is  the  Rayleigh 
probability  density  function  of  H(t)  (Appendix  A)  defined  by  Chalorabarti  (1987)  as 


p(H)  =  - — exp 
4mo 


-  Sm^  . 


(2.7) 


and  mo  is  the  variance  of  the  time  series.  If  the  following  substitution  is  made  (Sarpkaya 
and  Isaacson  1981), 

H„„=  5.09l/^ 

Eq.  (2.6)  may  be  solved  numerically  to  obtain  E[q;']  =0.001346. 
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CHAPTER  m:  Resolution  of  Incident  and  Reflected  Wave  Time  Series 


In  order  to  determine  experimental  wave  heights  and  to  calculate  envelope 
exceedance  coefficients,  the  incident  time  series  must  be  known.  Because  of  reflection 
from  models  in  a  2-dimensional  wave  channel,  a  method  to  resolve  the  incident  and 
reflected  random  wave  time  series  from  wave  gauges  must  be  used.  Bad  data  points  in 
wave  records  caused  by  wave  gauge  malfunction  are  smoothed  using  the  FORTRAN 
program  CLEANREC  that  employs  a  smoothing  spline  package  developed  by  H.  J. 
Woltring  (1986).  A  Fast  Fourier  Transform  (FFT)  analysis  is  used  to  compute  the  raw 
wave  spectrum  from  which  the  incident  and  reflected  random  wave  time  series  may  be 
determined.  Goda  and  Suzuki  (1976)  modified  a  method  developed  by  Thornton  and 
Calhoun  (1972)  to  resolve  the  incident  and  reflected  wave  spectra  from  wave  gauges 
spatially  separated  by  a  distance  Af  (see  Fig.  3.1). 


L( 


Figure  3.1  -  Definition  Sketch  for  two  wave  gauges 
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Kimura  (1985)  and  Fassardi  (1993)  extended  the  Goda  and  Suzuki  algorithm  to 


include  wave  phases  so  that  incident  and  reflected  time  series  could  be  resolved. 

Goda  and  Suzuki  (1976)  analyzed  two  simultaneous  wave  records  recorded  by 
closely  spaced  wave  gauges  aligned  in  the  direction  of  wave  propagation  by  an  FFT.  If 
the  incident  tj  and  reflected  f  random  wave  time  series  at  the  i*  wave  gauge  at  frequency 
CO  are  given  by 

T/j=  a  cos(kXj-  cot+  £)  (3- la) 

^.=  b  cos(kXj+  cot+  jS)  (3.1b) 

then  the  wave  profile  at  the  i*  wave  gauge  may  be  expressed  as 

(7?i+  ^'i)(x=x^=  A;coscot+  B;Sincot  (3.2) 

where 


Aj=  a  cos\/'j+  b  cosi/'j^  (3.3a) 

Bj=  a  simZ-j-  b  sin\j/^  (3.3b) 

A2=  a  cos(kAf  +  i/'j)+  b  cos(kAf  +  (3.3c) 

B2=  a  sin(kAf  +  b  sin(kAf  +  (3.3d) 

where  the  spatial  phases  of  the  incident  and  reflected  waves,  respectively,  are  defined  as 

i/',=  kXj+  £  (3.4a) 

kXj+  jS  (3.4b) 

X2=  Xj+  Af  (3.4c) 


and  Xi  and  X2  are  the  wave  gauge  positions  as  shown  in  Figure  3.1.  The  coefficients  A; 
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and  Bj  may  be  related  to  the  complex -valued  FFT  coefficient  X  by  (Dean  and  Dalrymple 
1984) 


where  j-  i/-  1  .  Goda  and  Suzuki  (1976)  solve  Eqs.  (3.3)  simultaneously  for  estimates 

of  the  amplitudes  a  and  b  listed  in  Table  3.1. 

For  irregular  waves,  the  time  series  are  the  superposition  of  many  wave 
components  given  by  Eqs.  (3.1).  Because  Eqs.  (3.3)  apply  to  each  Fourier  component  in 
an  irregular  wave  record,  linear  superposition  may  be  used  to  resolve  time  series  from  a 
spectrum  of  Fourier  components  (Goda  1985).  Goda  and  Suzuki  (1976)  observed  that 
spectral  estimates  diverged  (toward  infinity)  near  frequencies  satisfying  the  condition 
kAf  =n7r  for  n=0,l,2,...,  because  the  term  [  sin  kAf  |  in  the  denominator  of  the 
equations  for  the  amplitudes  a  and  b  in  Table  3.1  becomes  small  and  errors  from  noise 
are  amplified.  Therefore,  the  wave  gauge  spacing  Af  determined  the  frequency  limits  of 
a  band  pass  filter  from  which  the  waves  could  be  separated  into  incident  and  reflected 
time  series.  The  wave  amplitudes  a  and  b  could  be  resolved  effectively  for  the  interval 
0.  l7r<kAf  <0.97r.  Goda  and  Suzuki  (1976)  recommended  the  following  effective  band 
pass  limits  for  experimental  conditions: 

<Af<0.45L.  (3.6a) 

where  and  denote  the  wavelengths  corresponding  to  the  lower  (f^,^)  and  upper 
(fm.x)  frequency  limits,  respectively,  of  the  band  pass  filter.  Goda  and  Suzuki  (1976)  note 
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that  the  effective  range  may  be  taken  slightly  wider  than  given  in  Eq.  (3.6a)  such  that 

0-03L^^Af  <0.45L^  (3-6b) 

Goda  and  Suzuki  (1976)  also  recommended  that  the  wave  gauges  be  located  at  least  one 
wavelength  away  from  both  the  breakwater  toe  and  the  wave  generator. 

Kimura  (1985)  extended  the  Goda  and  Suzuki  algorithm  to  include  situations 
where  each  incident  wave  component  in  the  Fourier  spectrum  results  in  a  different 
reflection  coefficient  from  the  structure.  Kimura  (1985)  defined  incident  and  reflected 
wave  profiles  by  Eqs.  (3.1),  and  the  spatial  phases  by  Eqs.  (3.4a,b).  Eqs.  (3.1)  are  similar 
to  Kimura’s  Eqs.  (9  &  12)  on  pages  62-63  (1985)  if  the  sign  of  the  temporal  term  is 
changed  from  minus  to  plus  in  his  Eq.  (12).  Kimura  (1985)  also  solved  for  the  spatial 
wave  phases  defined  in  Eqs.  (3.4a,b).  Values  for  a,  b,  and  may  be  calculated  using 
the  equations  listed  in  Table  3.1  for  Kimura.  To  improve  the  resolution  of  incident  and 
reflected  waves,  Kimura  recommended  a  slightly  more  conservative  data  window  of 

0.15L  <Af<0.35L.  0-7) 

max  mm 

Fassardi  (1993)  followed  the  approach  of  Kimura  (1985)  to  extend  the  Goda  and 
Suzuki  algorithm  (1976).  However,  he  defined  the  incident  and  reflected  ^time  series 
at  the  i'^  wave  gauge  as 

?7j=  a  cos(kXj-  ut-  £)  (3.8a) 

b  cos(kx.+  wt+  iS)  (3.8b) 

If  the  spatial  phase  of  the  incident  and  reflected  waves,  respectively,  are  defined  as 

i/'j=  kx,-  £  (3-9a) 
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(3.9b) 


’/'r=  kx,"-  ^ 

rather  than  by  Eqs.  (3.4a,b),  then  the  amplitudes  and  spatial  phases  are  given  by  the 
equations  listed  in  Table  3.1  for  Fassardi  (1993).  The  only  difference  between  the 
Fassardi  and  Kimura  algorithm  is  the  sign  of  e  in  the  definition  of  the  incident  time 
series  phase.  This  sign  change  results  in  corresponding  sign  changes  in  the  Fourier 
coefficients  that  are  now  given  by 


A|=  a  cosi/',+  b  cos)/'^ 

(3.10a) 

Bj=  a  sim/'j+  b  sim/'^ 

(3.10b) 

A,=  a  cos(-kA£+  >/',)+  b  cos(kA£  + 

(3.10c) 

B,=  a  sin(-kA£+  b  sin(kA£  + 

(3.10d) 

These  four  equations  are  then  solved  simultaneously  to  obtain  the  equations  for  a,  b,  \pi, 
and  i//R  listed  in  Table  3. 1  for  Fassardi  (1993). 

In  the  experiments  documented  by  Medina  et  al.  (1994),  the  Fassardi  algorithm 
was  used  to  resolve  the  incident  and  reflected  time  series.  Three  sonic  wave  gauges 
aligned  in  the  direction  of  wave  propagation  were  centered  10m  from  the  toe  of  the 
breakwater  and  spatially  separated  by  Af  =  1.22m  according  to  Figure  3.2. 

Clean  time  series  from  each  gauge  were  used  in  the  FORTRAN  program 
RECANSYG  to  resolve  the  incident  and  reflected  time  series.  RECANSYG  outputs  the 
raw  time  series  and  spectra  resolved  from  the  wave  gauge  system  in  Figure  3.2.  The 
program  computes  the  time  series  mean  and  variance.  As  w'ell,  it  computes  the 
reflection  coefficient,  significant  wave  height,  and  mean  period  from  the  moments.  The 
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Table  3.1 -Summary  of  incident  and  reflected  wave  amplitudes  and  phases 
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incident  and  reflected  time  series  were  resolved  using  Gauges  1  and  2  and  Gauges  2  and 
3  so  that  the  algorithm  could  be  verified  by  duplication.  A  sample  run  is  analyzed.  As 
explained  in  APPENDIX  B  -  Simulation  of  Sea  States,  Run  E1P1L7  represents  a  sea 
state  simulated  with  the  following  parameters; 

El  =  Envelope  1  ;  7  =  10,  a;=2.52 
PI  =  Phase  1  ;  (^=27r/3 
L7  =  Height  7  ;  =0.723 


Figure  3.2 -Sonic  Wave  Gauge  Locations  (Fassardi  1993) 


For  a  small  data  sample  of  ten  seconds,  Figure  3.3  shows  how  the  clean  time 
series  are  phase  shifted  at  each  wave  gauge  station.  Figure  3.4  shows  two  minutes  of  the 
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clean  time  series  at  Gauges  1  and  2  and  the  resulting  incident  and  reflected  time  series 
resolved  by  the  Fassardi  (1993)  algorithm.  The  time  series  at  Gauge  1  is  identical  to 
that  at  Gauge  2,  except  for  the  phase  shift.  The  incident  and  reflected  time  series 
resolved  from  Gauges  1  and  2  are  plotted  on  the  same  scale  to  show  the  relative 
magnitude  of  the  incident  and  reflected  wave  amplitudes.  The  magnitude  of  the  wave 
amplitudes  of  the  reflected  time  series  are  about  25%  of  the  incident  amplitudes.  Figure 
3.5  shows  the  composite  spectrum  at  Gauge  2  and  the  incident  and  reflected  spectra 
resolved  from  Gauges  1  and  2  and  Gauges  2  and  3,  respectively.  The  spectra  in  Figure 
3.5  have  been  smoothed  using  the  non-statistical  technique  Box  Car  Averaging.  The 
resolved  spectra  show  that  the  energy  in  the  reflected  spectrum  is  about  5-10%  of  that  in 
the  incident  spectrum.  Breakwater  reflection  of  about  25%  observed  in  the  time  series 
results  in  the  reflected  spectral  density  being  5-10%  of  the  incident  spectral  density,  for 
each  frequency  component.  The  low  frequency  spike  in  the  spectral  density  at 
approximately  0.04  Hz  is  attributed  to  a  seiche  in  the  wave  channel  caused  by  the  start¬ 
up  of  the  wave  generator.  The  resolved  spectra  in  Figure  3.5  are  similar,  indicating  that 
the  Fassardi  (1993)  algorithm  produces  the  same  spectra  from  different  wave  gauge 
pairs.  At  the  truncation  limits  of  0.2140  Hz  and  0.7644  Hz,  the  spectral  density  of  the 
incident  and  reflected  waves  begins  to  diverge  because  the  term  |  sin  kAf  |  in  the 
denominator  of  the  equation  for  the  amplitudes  a  and  b  in  the  Fassardi  algorithm 
becomes  small.  Further  examples  of  the  resolution  of  incident  and  reflected  time  series 
and  spectra  from  the  output  of  RECANSYG  for  the  four  envelopes  discussed  in 
APPENDIX  B  are  included  in  APPENDIX  C  -  Time  Series  and  Spectra. 
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Time  Series  -  Run  E1P1L7 


-  -  GAUGE  1  —  GAUGE  2  -  -  -  GAUGE  3 


Figure  3.3  -  Phase  lag  between  time  series  at  three  wave  gauges 
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Figure  3.4  -  Time  series  resolved  from  Gauge  1  and  2 
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Figure  3.5  -  Composite  and  Resolved  Spectra 
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SUMMARY  AND  CONCLUSIONS 


Multiple  failures  of  existing  rubble  mound  breakwaters  designed  with  the  Shore 
Protection  Manual’s  (1984)  methodology  have  indicated  a  need  to  improve  the  design 
procedures  for  rubble  mound  structures.  Many  structural  and  environmental 
characteristics  such  as  armor  grading,  structure  permeability,  wave  period,  storm 
duration,  spectral  shape,  and  wave  groupiness  have  secondary  effects  on  breakwater 
stability.  Irregular  waves  generated  solely  from  an  energy  spectrum  with  random  phases 
do  not  accurately  model  real  wave  trains,  and  the  succession  of  wave  heights  or 
groupiness  is  an  important  factor  in  the  design  of  rubble  mound  breakwaters.  The 
magnitude  of  the  energy  flux  given  by  the  envelope  exceedance  coefficient  a  is  one  way 
of  quantifying  the  groupiness  of  a  wave  time  series.  Intermediate  scale  wave  channel 
studies  at  Oregon  State  University  documented  by  Medina  et  al.  (1990),  Fassardi  (1993), 
and  Medina  et  al.  (1994)  have  demonstrated  that  the  envelope  exceedance  coefficient  a 
may  be  used  to  correlate  wave  groupiness  with  breakwater  damage.  No  significant 
correlation  was  found  between  the  peak  enhancement  factor  y  and  rubble  mound 
breakwater  stability.  Furthermore,  the  constant  phase  shift  applied  to  the  time  series 
for  each  wave  envelope  A(x,t)  did  not  effect  the  armor  damage. 

To  calculate  the  envelope  exceedance  coefficient  a  from  the  incident  wave  time 
series,  either  the  Kimura  (1985)  or  the  Fassardi  (1993)  algorithm  may  be  used  to  resolve 
the  incident  and  reflected  waves  from  closely  spaced  wave  gauges.  The  Fassardi  (1993) 
algorithm,  used  by  Medina  et  al.  (1990)  in  their  experiments  at  Oregon  State  University 
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is  different  from  the  Kimura  (1985)  algorithm  by  the  sign  of  the  random  phase  angle  of 
the  incident  time  series  (Eq.  3.1a  &  3.8a).  This  sign  change  does  not  have  any  effect  on 
the  value  of  a  or  the  interpretation  of  the  results.  For  the  frequency  range  given  by 
Table  B.l,the  Fassardi  (1993)  algorithm  effectively  resolves  incident  and  reflected  wave 
time  series  from  closely  spaced  wave  gauges  aligned  in  the  direction  of  wave 
propagation.  Some  divergence  of  the  incident  and  reflected  wave  spectra  is  observed 
near  the  cut  off  frequencies  of  the  truncated  spectra. 

Since  the  wave  resolution  algorithms  presented  only  resolve  wave  components  in 
the  frequency  range  determined  by  the  wave  gauge  spacing,  the  spectra  used  to  simulate 
the  time  series  for  the  laboratory  experiments  were  truncated  accordingly.  The  variance 
of  the  full  spectra  was  preserved  for  the  truncated  spectra  by  considering  the  squared 
significant  wave  height  ratio  given  by  Eq.  (B.6).  For  the  physical  simulations  performed 
by  Medina  et  al.  (1990),  phase  spectra  resulting  in  the  highest  and  lowest  values  for  the 
envelope  exceedance  coefficient  a  were  chosen  for  each  amplitude  spectra.  Significant 
wave  heights  of  consecutive  runs  were  increased  in  discrete  increments  such  that  the 
stability  numbers  for  the  two  armor  rock  sizes  were  equal  in  consecutive  runs. 
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APPENDIX  A:  Rayleigh  Distribution 


The  Rayleigh  probability  density  function  (pdf)  for  the  wave  amplitudes  a  (  = 
H/2)  is  given  by 


P(a)  =  —exp 
"r 


1 

2 


a 


U(a)  ;  o;j^>0 


(A.  la) 


where  =  Rayleigh  parameter  (Hoffman  and  Karst  1975)  and  U(a)  =  Heaviside  step 
function.  The  Rayleigh  cumulative  distribution  function  (cdf)  for  wave  amplitudes  is 


P(a)  =  1  -  exp 


1 

2 


a 

Oir, 


U(a) 


The  following  4  forms  of  (A.  la)  may  be  found: 


(A.  lb) 


7ra 


p(a)  =  ^exp 


2H 

P(H)  =  -^exp 

H 


p(a)  =  — exp 
m. 


.TT^  H 

p(H)  =  - — exp 
4m 

o 


TT 

T 


a 


2 

H 

H 

rms 

1 

2 

J. 

2 


U(a)  ; 


U(H)  ;  H  >0 

^  '  '  rms 


U(a)  ;  m  >0 


H 


2\/m 


U(H) 


m  >0 

o 


(A.  2a) 


(A.2b) 


(A.  2c) 


(A.  2d) 


where  /i,  =  average  of  the  amplitude  a;  H^,  =  root-mean-square  wave  height;  and  m,,  = 
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variance  of  the  time  series  for  the  water  surface  elevation  7j(t).  When  the  random 
variable  a  in  (A. 2)  is  scaled  by  the  standard  deviation  of  the  time  series  /in^,  then  the 
following  change  of  variables  gives; 


p(?)  =  |exp 


(A.2e) 


Eqs.  (A. 2)  may  be  derived  from  (A.  la)  by  solving  for  the  generic  Rayleigh  parameter 
The  percent  or  fraction  of  wave  amplitudes  greater  than  wave  amplitude 
a„  (=  Hn/2),  may  be  computed  from  (A. la)  by 


CO 


and  the  natural  logarithm  of  (A. 3)  gives 
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=  \/-21n(n)  (A.4) 

"r 


where  0  <  n  <  1 .  The  average  wave  amplitude  greater  than  a^  may  be  determined  from 


p(a)da  =  J  ap(a)da 
that  may  be  integrated  by  parts  to  obtain 


2 

CD 

a  (n)  =  na  i/-21n(n)  +  [exp 

_  1 

a 

J 

2 

.. 

“r 

by  (A. 3)  and  (A.4).  Finally,  integration  gives 


IL  =  ^/-21n(n)  +  1 
n  n 


T 

T 


1  -erf 


(/-ln(n)) 


(A.  5a) 


=  v/-21n(n) 


X 


erfc 


(v/-ln(n)  ) 


(A.5b) 


where  erf(  • )  =  the  error  function  and  erfc(  • )  =  complementary  error  function 

(Kreyszig  1983).  In  order  to  relate  the  Rayleigh  parameter  q:r  to  the  average  /x.,  set  n  = 
1  in  Eq.  (A. 5)  and 


and  substitution  into  (A.  la)  gives  (A. 2a). 


(A.6) 
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Alternatively,  (A. 5)  may  be  derived  directly  from  special  functions  according  to 

oo  00 

i  I  p(a)da  =  |  ap(a)da 


q:„ 


-  -  v^r 

2 


3  a ' 

_  0  “ 
—  j'J ) 


2  2.^, 


where  the  generalized  incomplete  Gamma  function  r(a,Zo,Zi)  (Wolfram  1991)  is  defined 


as 


r(a,z„.z,)  -  r{a.2<,)  -  r(a,z,) 

and  the  incomplete  Gamma  function  r(a,z)  (Wolfram  1991)  is  defined  as 

oa 

r(a,z)  =  r(a)  -  G(a,z)  =  |  t*'‘exp( -t)  dt 


G(a,z)  = 


t  *'‘exp(-t)  dt 


Substituting  for  r(a,Zo,Zi) 
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3  a' 

U  n 


2  2-^ 


On 
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=  r 


r  3 
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v/T exp(  -t)  dt 


i/r exp(  -t)  dt 


=  n\/-ln(n)  +  1  -erf(\/ -In (n)  ) 

=  n\/-ln(n)  +  2^erfc(\/-ln(n)  ) 
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giving  (A. 5). 

To  relate  the  Rayleigh  parameter  o;r  to  make  the  following  change  of 
variables  in  (A.  la): 


y  =  a^  ]  dy  =  2ada  -  2\fyd‘2. 


and 


p(a)da  =  p(y)dy 


giving 


P(y) 
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exp 


_  y 


U(y) 


(A.7) 


The  average  of  a^  may  now  be  computed  from 


a^  = 


'mis 


I  yp(y)dy 


CO 


(A.  8) 


and  the  Rayleigh  parameter  o!r  may  be  defined  by 
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(A.9) 


« 


H. 


«R  = 


v/2  2^2 


Substituting  (A. 6)  and  (A.9)  into  (A. 5b)  gives 


H 


Q 


\/-ln(n) 


y-81n(n) 


4- 


efrc (\/-ln(n)  ) 
n 


The  significant  wave  height  or  average  of  the  highest  1/3  wave  heights  I 
computed  from  (A.5b)  for  n  =  1/3  (dropping  the  overbar  average  notation) 

^  ^  =  2.00215 

2  Or  2  Or 

2.00215 


Substituting  (A. lib)  into  (A. 5b)  gives 


(A.  10a) 

(A. 10b) 

(A.  10c) 

1/3  may  be 
according  to 

(A.  11a) 

(A. 11b) 
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H 


=  (2.00215)-^ 


(A.  12) 
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1  erfc(v/-ln(n)  ) 


To  find  the  mode  (or  most  probable)  value  of  a^^de  (=  H„odc/2),  the  maximum  value  of 
(A.  la)  occurs  when 


dp(a)  ^  Q 
da 


^  ^mode 


(A.  13) 
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(A.15a,b) 


H.pws) 


(2.00215)-’ 


(A.15c,d) 


The  mean  or  median  of  a  occurs  when  P(a)  =  0.5  in  (A.  lb).  The  natural  logarithm  of 
(A.  lb)  for  P(a„cd)  =  0.5  is 

=  s/2\n{2)  (A.  16) 

“r 
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giving 


^mcd 


41n(2) 


(A.17a,b) 


=  v/81n(2) 


H(H,,3) 


v/21n(2) 

2.0015 


(A.17c,d) 


These  results  are  summarized  in  Table  A.l 
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4.06198 

1.79992 

5.09094 
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5.32423 

2.35924 

6.67293 

1.66644 

Table  A.l -Rayleigh  distributed  wave  height  relationships 
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APPENDIX  B:  Simulation  of  Sea  States 


In  order  to  correctly  describe  wave  groups  in  laboratory  studies,  the  incident  wave 
train  must  be  known.  An  algorithm  developed  at  Oregon  State  University  by  Fassardi 
(1993)  was  used  to  resolve  the  incident  and  reflected  wave  trains  in  a  finite  bandwidth 
determined  by  the  wave  gauge  spacing  Af  (Table  3.1).  One  broad  (7  =  !)  and  one 
narrow  (7  =  10)  truncated  Goda-JONSWAP  spectrum  (Goda  1985)  were  simulated 
according  to 


s,(f)= 


I  R  s,lninc 


f*  f 

^exp[-1.25(J-)- 


It 


cxp[- 


(B.l) 


0.0624 

0.230+  0.03367-  0. 185(1.9+  7)-’ 


cr= 


0.07  ;  f<f 
0.09  :  f>f 

P 


(B.2) 

(B.3) 


for  f,n^  <  f  <  fn,jx,  where  fp  is  the  peak  frequency,  Hr  is  the  squared  significant  wave  height 
ratio  defined  in  Eq.  (B.6),  and  is  the  significant  wave  height  of  the  truncated 

spectrum.  Eq.  (B.l)  is  different  than  Eq.  (30)  in  Fassardi  (1993)  in  that  Fassardi  (1993) 
does  not  document  the  use  of  the  constant  Ci  defined  by  Eq.  (B.2).  The  spectra  were 
truncated  according  to 


f  .  =0.7f 

imn  p 

f  =2.5f 

max  p 


(B.4) 
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To  define  the  truncated  spectra  and  preserve  the  variance  of  the  full  spectra,  Eq. 

(B.  1)  includes  the  dimensionless  coefficient  Hr.  If  the  significant  wave  height  is  given  in 
terms  of  the  variance  (Table  A.l)  by 

H^=4.0043/m^ 

then  the  squared  significant  wave  height  ratio  Hr  may  be  defined  as 

u  _  ^s.fuu  (B.6) 

R  H^ 

where  H,  f^u  is  the  significant  wave  height  of  the  full  spectra.  For  peak  frequencies 
between  0.1  and  0.5  Hz,  the  average  dimensionless  variance  preserving  coefficients  for 
7  =  1  and  7  =  10  are  Hr=  1.0375  and  Hr=  1.0133  respectively. 

The  mean  frequency  of  the  truncated  spectra  defined  by  the  first  and  zeroeth 

spectral  moments 

(B.7) 

peak  frequency  of  the  truncated  spectra 

(B.8) 

For  peak  frequencies  between  0.1  and  0.5  Hz,  the  average  frequency  ratios  were 


was  held  constant  at  1=0.3333  Hz  such  that  the 
could  be  computed  from  the  frequency  ratio 


f 

f  =-l 
^  T 
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4=0.81 15  and  4=0.9173  for  7  =  1  and  7=10,  respectively.  Accordingly,  the  peak 
frequencies  used  to  simulate  the  truncated  spectra  from  Eq.  (B.l)  were  fp=0.2705Hz 
and  fp=0. 3058Hz  for  7  =  1  and  7  =  10,  respectively. 

The  wave  gauge  spacing  defined  by  Figure  3.1  was  1.22m,  and  the  water  depth 
was  3.05m  (Fassardi  1993).  If  the  band  pass  filter  limits  proposed  by  Goda  and  Suzuki 
(1976)  given  by  Eq.  (3.6a)  are  used,  and  the  linear  wave  theory  dispersion  relationship  is 
given  by 

the  minimum  and  maximum  resolvable  frequencies  for  Af  =  1.22m  are  4^=0. 1299  Hz 
and  4,x=0.7589  Hz.  For  cases  where  7  =  1,  the  minimum  and  maximum  frequencies  of 
the  truncated  spectra  are  4jn=0.1894Hz  and  43,^=0.6762  Hz.  For  7  =  10,  the  minimum 
and  maximum  frequencies  of  the  truncated  spectra  are  4^=0. 2140Hz  and  4„=0.7644 
Hz.  These  frequency  limits  for  the  truncated  spectra  are  reasonably  close  to  those 
resolvable  by  wave  gauges  spatially  separated  by  Af  =  1.22m. 

Different  phase  spectra  produce  time  series  with  different  wave  grouping 
characteristics.  Accordingly,  time  series  were  synthesized  using  two  different  phase 
spectra  for  each  of  the  truncated  Goda-JONSWAP  spectra  given  by  Eq.  (B.l)  with 

7  =  l,10and  ?=  0.3333  Hz.  The  phase  spectra  were  chosen  from  100  random 

Deterministic  Spectral  Amplitude  (DSA)  simulations  for  each  amplitude  spectra.  Values 
of  a  were  calculated  from  each  DSA  simulated  time  series  according  to  Eq.  (2.5).  For 


Ilf  ,  2'?rh 

g — tanh - 

°  L  L 


(B.9) 
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each  7,  the  phase  spectra  that  produced  the  wave  envelope  A(x,t)  with  the  highest  and 
lowest  wave  height  variability  characterized  by  a  were  chosen  for  the  physical 
simulations.  The  parameters  used  for  wave  simulation  are  summarized  in  Table  B.l. 


7 

Cl 

Hk 

fp 

resolvable  for 
Af  =  1.22m 
from  Eq.  (3.6a) 

f  . 

•^min 

f 

•^max 

Low  a 
from 
DSA 

High  a 
from 
DSA 

f  . 

f 

•^max 

1 

0.3123 

1.0375 

0.8115 

0.2705 

0.1299 

0.7589 

0.1894 

0.6762 

0.51 

2.18 

10 

0.1134 

1.0133 

0.9173 

0.3058 

0.1299 

0.7589 

0.2140 

0.7644 

0.23 

2.52 

Table  B.l  -  Summary  of  simulation  parameters  for  Af  =  1.22m  and  h=3.05m 

(frequencies  in  Hz) 


Realizations  with  the  same  envelope,  but  different  phase-shifted  wave  profiles  are 
given  by  (Medina  et  al.  1994) 


M 

’?i(^>0  =  X^a^cos[kx-(jt-(e+(^)] 

m  =  I 


(B.IO) 


where  0  is  a  constant  phase  shift  applied  to  all  M  wave  components.  Three  different 
phase  shifted  realizations  (<^>=0,2-71/3,  47r/3)  and  the  replicate  (4>=2'jr)  for  each  of  the 
four  wave  envelopes  were  tested  to  determine  if  realizations  shifted  by  a  constant  phase 
produce  different  levels  of  damage.  For  each  y,  A(x,t),  and  constant  phase  shift  4>,  a 
sequence  of  time  series  were  tested  with  increasing  significant  wave  heights. 
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The  breakwater  model  tested  by  Medina  et  al.  (1990)  was  divided  in  the  direction 
of  wave  propagation  into  two  sides  with  a  different  armor  layer  rock  size  (viz.,  Wl=  128.5 
N  and  Ws=99.1  N)  on  each  side.  The  monochromatic  wave  height  corresponding  to  the 

zero  damage  condition  (SPM  1984)  can  be  calculated  from 

1 

■j 

(B.ll) 


3 

WK, 

A-i 

cot6 

P. 

Pr 

where  W  is  the  median  value  of  the  mass  distribution  of  rocks  in  the  armor;  is  the 
stability  coefficient;  0  is  the  angle  of  the  breakwater  slope  measured  from  horizontal, 
and  Pr  and  p^.  are  the  weight  densities  of  the  rocks  and  water,  respectively.  Solving  Eq. 
(B.ll)  numerically  for  the  small  rock  size  (Ws=99.1  N)  yields  a  monochromatic  wave 
height  of  H=0.55m  for  the  zero  damage  condition.  Breakwater  model  and  armor  rock 
characteristics  are  given  in  Table  B.2. 


W,  (N) 

Ws  (N) 

Kd 

p,  (kN/m^) 

p„  (kN/m^) 

COt0 

128.5 

99.1 

4 

21 A 

9.8 

2 

Table  B.2  -  Summary  of  breakwater  model  characteristics  (Fassardi  1993) 


For  the  zero  damage  condition,  if  the  representative  design  wave  height  for  a 
random  sea  is  Hj/io,  then  the  design  significant  wave  height  for  the  small  rock  size  is 
determined  (Table  A.l)  by 
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H  0.55m  ^  for  W  =99.1  N 
'  1.27137  " 


(B.12) 


An  empirically  derived  number  which  controls  the  stability  of  the  armor  layer  (Medina  et 
al.  1990)  is  given  by 


where  Nj  is  the  stability  number,  and  is  the  design  significant  wave  height.  The 
significant  wave  heights  in  the  physical  simulations  were  increased  in  discrete  increments 
such  that  the  stability  numbers  given  by  Eq.  (B.13)  for  the  large  (WJ  and  small  (Wg) 
rocks  were  equal  in  consecutive  runs.  Seven  significant  wave  heights  were  tested  for 
each  of  the  sixteen  realizations.  The  significant  wave  heights  varied  from  the  design 
significant  wave  height  of  the  small  rock  size  to  the  maximum  wave  height,  avoiding 
breaking.  From  Eq.  (B.13),  the  significant  wave  height  of  the  k*  run  is  given  by 


W, 

H  =0.43  -J: 

s.Oc)  ^ 


k=  1,  2, 


(B.14) 


Thus,  twenty-eight  realizations  of  four  phases  and  seven  wave  heights  were 
simulated  for  each  of  four  target  envelopes  shown  in  Figure  B.l.  For  each  realization, 
runs  of  approximately  30  minutes  long  containing  N=2’^=32768  points  sampled  at 
At=0.06  seconds  were  simulated. 
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Figure  B.l  -  Simulation  of  Realizations  with  Wave  Height 


APPENDIX  C:  Wave  Time  Series  and  Spectra 


Section  Ease 

1.  Run  E1P1L7 . 38 

2.  Run  E2P1L7 . 40 

3.  Run  E3P1L7 . 45 

4.  Run  E4P1L7 . 50 
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RESOLVING  INCIDENT  AND  REFLECTED  TRAINS  FOR  RECORD:  elplIV 


Gauge  No.  :  I 

Raw  Time  Series  Mean  Value  :  0.009327  Variance  : 

Gauge  No.  :  2 

Raw  Time  Series  Mean  Value  :  0.010709  Variance  : 
Gauge  No.  :  3 

Raw  Time  Series  Mean  Value  :  0.006861  Variance  : 

RESOLVED  INCIDENT  AND  REFLECTED  WAVES 

Analysis  from  Gauges  2  and  3 

Incident  wave  Spectrum’s  Variance:  0.032890 
Reflected  wave  Spectrum’s  Variance:  0.001678 

Breakwater’s  reflection  coefficient  :  0.2259 
Incident  significant  wave  height  :  0.7261 
Reflected  significant  wave  height  :  0.1640 
Mean  Period  (TOl)  from  incident  spectrum  :  3.065 

Incident  wave  train,  Mean  and  Variance  from  TS(2,3) 

TS  Mean  :  0.000000  TS  Variance  :  0.032890 

Reflected  wave  train,  Mean  and  Variance  from  TS(2,3) 
TS  Mean  :  0.000000  TS  Variance  :  0.001678 

Analysis  from  Gauges  1  and  2 

Incident  wave  Spectrum’s  Variance:  0.033244 
Reflected  wave  Spectrum’s  Variance:  0.002225 

Breakwater’s  reflection  coefficient  :  0.2587 
Incident  significant  wave  height  :  0.7300 
Reflected  significant  wave  height  :  0.1889 
Mean  Period  (TOl)  from  incident  spectrum  :  3.074 

Incident  wave  train.  Mean  and  Variance  from  TS(1,2) 

TS  Mean  :  0.000000  TS  Variance  :  0.033244 

Reflected  wave  train.  Mean  and  Variance  from  TS(1,2) 
TS  Mean  :  0.000000  TS  Variance  :  0.002225 


0.040746 

0.038375 

0.036691 
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RESOLVING  INCIDENT  AND  REFLECTED  TRAINS  FOR  RECORD:  e2pll7 


Gauge  No.  :  1 

Raw  Time  Series  Mean  Value  :  0.009919  Variance  : 

Gauge  No.  :  2 

Raw  Time  Series  Mean  Value  :  0.003559  Variance  : 

Gauge  No.  :  3 

Raw  Time  Series  Mean  Value  :  -0.010094  Variance  : 

RESOLVED  INCIDENT  AND  REFLECTED  WAVES 

Analysis  from  Gauges  2  and  3 

Incident  wave  Spectrum’s  Variance;  0.031342 
Reflected  wave  Spectrum’s  Variance:  0.002484 

Breakwater’s  reflection  coefficient  :  0.2815 
Incident  significant  wave  height  :  0.7089 
Reflected  significant  wave  height  :  0.1996 
Mean  Period  (TOl)  from  incident  spectrum  ;  3.245 

Incident  wave  train,  Mean  and  Variance  from  TS(2,3) 

TS  Mean  ;  0.000000  TS  Variance  :  0.031342 

Reflected  wave  train,  Mean  and  Variance  from  TS(2,3) 
TS  Mean  :  0.000000  TS  Variance  :  0.002484 

Analysis  from  Gauges  1  and  2 

Incident  wave  Spectrum’s  Variance:  0.033780 
Reflected  wave  Spectrum’s  Variance:  0.002935 

Breakwater’s  reflection  coefficient  ;  0.2948 
Incident  significant  wave  height  :  0.7359 
Reflected  significant  wave  height  ;  0.2169 
Mean  Period  (TOl)  from  incident  spectrum  :  3.244 

Incident  wave  train.  Mean  and  Variance  from  TS(1,2) 

TS  Mean  :  0.000000  TS  Variance  ;  0.033780 

Reflected  wave  train.  Mean  and  Variance  from  TS(1,2) 
TS  Mean  :  0.000000  TS  Variance  :  0.002935 


0.039530 

0.041591 

0.041708 
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Resolved  from  G 


Tinie  Ser-ies  -  JR.un 
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Time  Series  -  Run  E2P1L7 


RESOLVING  INCIDENT  AND  REFLECTED  TRAINS  FOR  RECORD:  e3plI7 


Gauge  No.  :  1 

Raw  Time  Series  Mean  Value  :  0.011743  Variance  : 
Gauge  No.  :  2 

Raw  Time  Series  Mean  Value  :  0.002714  Variance  : 
Gauge  No.  :  3 

Raw  Time  Series  Mean  Value  :  -0.002450  Variance  : 

RESOLVED  INCIDENT  AND  REFLECTED  WAVES 

Analysis  from  Gauges  2  and  3 

Incident  wave  Spectrum’s  Variance:  0.035032 
Reflected  wave  Spectrum’s  Variance:  0.001793 

Breakwater’s  reflection  coefficient  :  0.2262 
Incident  significant  wave  height  :  0.7494 
Reflected  significant  wave  height  :  0.1695 
Mean  Period  (TOl)  from  incident  spectrum  :  3.056 

Incident  wave  train,.  Mean  and  Variance  from  TS(2,3) 

TS  Mean  :  0.000000  TS  Variance  :  0.035032 

Reflected  wave  train,  Mean  and  Variance  from  TS(2,3) 
TS  Mean  :  0.000000  TS  Variance  :  0.001793 

Analysis  from  Gauges  1  and  2 

Incident  wave  Spectrum’s  Variance:  0.034591 
Reflected  wave  Spectrum’s  Variance:  0.002552 

Breakwater’s  reflection  coefficient  :  0.2716 
Incident  significant  wave  height  :  0.7447 
Reflected  significant  wave  height  :  0.2023 
Mean  Period  (TOl)  from  incident  spectrum  ;  3.078 

Incident  wave  train.  Mean  and  Variance  from  TS(1,2) 

TS  Mean  :  0.000000  TS  Variance  :  0.034591 

Reflected  wave  train.  Mean  and  Variance  from  TS(1,2) 
TS  Mean  :  0.000000  TS  Variance  :  0.002552 


0.041173 

0.038172 

0.038333 
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Tin-»e  Series  —  I?.xm  E3DPJLI-<'7 
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RESOLVING  INCIDENT  AND  REFLECTED  TRAINS  FOR  RECORD:  e4pll7 


Gauge  No.  :  1 

Raw  Time  Series  Mean  Value  :  0.015713  Variance  : 
Gauge  No.  :  2 

Raw  Time  Series  Mean  Value  :  0.025821  Variance  ; 

Gauge  No.  :  3 

Raw  Time  Series  Mean  Value  :  0.029771  Variance  : 

RESOLVED  INCIDENT  AND  REFLECTED  WAVES 

Analysis  from  Gauges  2  and  3 

Incident  wave  Spectrum’s  Variance:  0.036691 
Reflected  wave  Spectrum’s  Variance:  0.002226 

Breakwater’s  reflection  coefficient  :  0.2463 
Incident  significant  wave  height  :  0.7670 
Reflected  significant  wave  height  :  0.1889 
Mean  Period  (TOl)  from  incident  spectrum  :  3.206 

Incident  wave  train,  Mean  and  Variance  from  TS(2,3) 

TS  Mean  :  0.000000  TS  Variance  :  0.036691 

Reflected  wave  train,  Mean  and  Variance  from  TS(2,3) 
TS  Mean  :  0.000000  TS  Variance  :  0.002226 

Analysis  from  Gauges  1  and  2 

Incident  wave  Spectrum’s  Variance:  0.036596 
Reflected  wave  Spectrum’s  Variance:  0.002472 

Breakwater’s  reflection  coefficient  :  0.2599 
Incident  significant  wave  height  :  0.7660 
Reflected  significant  wave  height  :  0.1991 
Mean  Period  (TOl)  from  incident  spectrum  :  3.212 

Incident  wave  train.  Mean  and  Variance  from  TS(1,2) 

TS  Mean  :  0.000000  TS  Variance  :  0.036595 

Reflected  wave  train.  Mean  and  Variance  from  TS(1,2) 
TS  Mean  :  0.000000  TS  Variance  :  0.002472 


0.043060 

0.041359 

0.044257 
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APPENDIX  E:  Notation 


A(x,t) 

Ai 


^Tnc4 


^modc 


a 


n 


^rms 

a 

Bi 

b 

Cl 

da 

dH 

de 

E(a') 
erf(  • ) 
erfc(  • ) 
f 
I 

f 

•^mAx(min) 


=  wave  envelope  function 
=  FFT  coefficient  at  i'^  wave  gauge 
=  the  median  of  a 
=  most  probable  value  of  a 

=  the  average  wave  amplitude  greater  than  a„ 

=  root  mean  square  wave  amplitude 
=  incident  wave  amplitude 
=  FFT  coefficient  at  i*^  wave  gauge 
=  reflected  wave  amplitude 

=  Goda-JONSWAP  spectrum  parameter  defined  in  Fq.  (B.2) 
=  wave  amplitude  interval 
=  wave  height  interval 
=  interval  of  ^  defined  in  Fq.  (A.2e) 

=  expected  value  of  a' 

=  error  function  of  (  • ) 

=  complementary  error  function  of  (  • ) 

=  wave  frequency 
=  mean  frequency  of  spectrum 
=  maximum  (minimum)  cut-off  frequency 
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G(a,z) 


H(x,t) 


H(nAt) 


=  peak  frequency  of  spectrum 
=  frequency  ratio  given  by  Eq.  (B.8) 

=  gamma  function  from  APP.  A 
=  gravitational  constant 
=  wave  height 

=  wave  height  function  at  a  fixed  location 
=  wave  height  function 
=  discrete  wave  height  function  at  x=0 
=  the  median  of  H 
=  the  most  probable  value  of  H 

=  squared  significant  wave  height  ratio  given  by  Eq.  (B.5) 

=  root  mean  square  wave  height 
=  significant  wave  height  =  H,/3 
=  significant  wave  height  of  full  spectrum 
=  significant  wave  height  of  the  k'*’  run 
=  significant  wave  height  of  truncated  spectrum 
=  average  of  the  wave  heights 

=  average  of  the  highest  one-third  of  the  wave  heights 
=  average  of  the  highest  one-tenth  of  the  wave  heights 
=  average  of  the  highest  one-hundredth  of  the  wave  heights 
=  water  depth 
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Ko 

k  (=27r/L) 
L 

Lmax(miii) 

ITlo 

mi 

n 

N 

N3 

P(-) 

P(  •) 

s,(f) 

T 


=  armor  layer  stability  coefficient 
=  wave  number 
=  wave  length 

=  maximum  (minimum)  wave  length 
=  variance  of  the  time  series,  zero-th  spectral  moment 
=  first  spectral  moment 

=  percent  or  fraction  of  wave  amplitudes  greater  than  a„ 
=  number  of  data  points 
=  stability  number  defined  in  Eq.  (B.12) 

=  Rayleigh  cumulative  distribution  function  of  (  • ) 

=  Rayleigh  probability  density  function  of  ( • ) 

=  Goda-JONSWAP  spectrum  defined  by  Eq.  (B.l) 

=  wave  period 


t 


=  time 


U(  • )  =  Heaviside  step  function  of  ( • ) 

W  =  the  median  value  of  the  mass  distribution  of  rocks  in  the  armor 

=  large  armor  layer  rock  size 
Ws  =  small  armor  layer  rock  size 

X  =  complex  valued  FFT  coefficient 

Xj  =  position  of  i'^  wave  gauge  from  wave  maker 

y  =  change  of  variable  from  Eq.  (A. 7) 

a  =  envelope  exceedance  coefficient 
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t  • 


“r 

IS 

r(a,z) 

r(a,Zo,Zi) 

T 

AH„ 

A£ 

At 


E 

Vi 

e 


Ma 

Pr 

Pw 

O’ 


=  magnitude  of  the  variation  of  wave  energy  above  Hj/jo 
=  Rayleigh  parameter 
=  reflected  wave  phase 
=  incomplete  gamma  function  from  APR.  A 
=  generalized  incomplete  gamma  function  from  APP.  A 
=  peak  enhancement  factor 

=  normalized  measure  of  the  variation  of  wave  height 

=  wave  gauge  spacing 

=  sampling  time  interval 

=  incident  wave  phase 

=  reflected  time  series  at  i*  wave  gauge 

=  sea  surface  elevation 

=  incident  time  series  at  i*  wave  gauge 

=  angle  of  the  breakwater  slope  measured  from  horizontal 

=  average  of  the  amplitude  a 

=  change  of  variable  from  Eq.  (A.2e) 

=  weight  density  of  rock 
=  weight  density  of  water 

=  Goda-JONSWAP  spectrum  parameter  defined  in  Eq.  (B.3) 


=  constant  phase  shift  applied  to  each  realization 
=  incident  (reflected)  spatial  wave  phase 
CO  (=27r/T)  =  circular  wave  frequency 
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